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In a number of current experiments in the field of spin-caloritronics a temperature gradient across a nanos¬ 
tructured interface is applied and spin-dependent transport phenomena are observed. However, a lack in the 
interpretation and knowledge let it unclear how the temperature drop across a magnetic nanostructured inter¬ 
face looks like where both phonons and electrons may contribute to thermal transport. We answer this question 
for the case of a magnetic tunnel junction (MTJ) where the tunneling magneto Seebeck effect occurs. Never¬ 
theless, our results can be extended to other nanostructured interfaces as well. Using an ab initio method we 
explicitly calculate phonon and electron thermal conductance across Fe/MgO/Fe-MTJs by using Green's func¬ 
tion method. Further, by estimating the electron-phonon interaction in the Fe leads we are able to calculate the 
electron and phonon temperature profile across the Fe/MgO/Fe-MTJ. Our results show that there is an electron- 
phonon temperature imbalance at the Fe-MgO interfaces. In consequence, a revision of the interpretation of 
current experimental measurements might be necessary. 

PACS numbers: 


In the 1990s, Johson and Silsbee |T1] have pioneered the the¬ 
oretical investigation of thermal effects in nano-magnetic sys¬ 
tems. However, it took until 2008 for spin-caloritronics Si 
phenomenon to emerge after the discovery of the spin Seebeck 


effect in NiFe films by Uchida et. al |4f 5]. Since then in¬ 
teresting spin-caloritronic transport phenomena are observed 
including thermal spin injection |0] and magneto-Peltier ef¬ 
fect (7}-[9]. Magneto-thermoelectric phenomena, e.g. the tun¬ 
neling magneto-Seebeck (TMS) effect, which occurs in mag¬ 
netic tunnel junctions (MTJs) and is the dependence of the 
charge Seebeck coefficient on the magnetic alignments, is 


demonstrated both in experiments |1( 


and first-principles 


calculations Efl. Also thermal spin transfer torque has 
been theoretically predicted | F5 ,:ll and experiments are on 
the way & This arise the possibility of MgO-MTJs as a 
potential key structure for spin-caloritronic applications. 


All these phenomena have in common that there is a tem¬ 
perature gradient across an interface or across a very thin 
layer. However, at present the understanding of the temper¬ 
ature drop and even how heat is transported across a magnetic 
interface is incomplete or unknown. In this letter we focus 
on MgO-based-MTJs, but our results can at least qualitatively 
applied to other systems as well. 

To get the correct temperature drop across the barrier in a 
MTJ the appropriate value of the thermal conductance n AITJ 
is important. Considering the fact that heat is carried by 
phonons as well as electrons it is not clear if n MTJ is just the 
sum of both contributions as in the clasical limit or if we need 
a new definition. In the analysis of the experiments the tem¬ 
perature across the MgO-MTJs has been simulated by COM- 
soleo where n MTJ is typically taken from the thermal 
conductivity of bulk-MgO or MgO thin film or in between, 
where the variation of this value is one order of magnitude. 


In order to check if this approximation is correct one should 
get a rough idea about the phonon transport length scale in 
Fe/MgO/Fe-MTJs. For this purpose, the phonon mean free 


path A of bulk MgO can be estimated by 1811 


k p = 1/3 CvA (1) 

where k p is the thermal conductivity, C is the volumetric spe¬ 
cific heat, and v is the average phonon group velocity. By tak¬ 
ing experimental values of these quantities for MgO, which 
are listed in Table Q] the phonon mean free path A for bulk 
MgO is around 6.4 nm. Consequently, A is larger than the 
typical MgO barrier thickness d (between 0.6~2 nm) used in 
MTJs. This means it is not appropriate to take for the thermal 
conductivity of the junction the bulk value of MgO or the thin 
film value because these thin films have a thickness of several 
hundreds of nanometer 0. Eventually, the vibrational prop¬ 
erties of the interface region and the phonon transport across 
the whole Fe/MgO/Fe junction has to be calculated explicitly. 

Further, thinking about a MTJ it is not obvious what is the 
role of electrons versus phonons for the thermal transport. In 
the metallic ferromagnetic leads the thermal transport is dom¬ 
inated by electrons whereas in MgO as an insulator the ther¬ 
mal transport is determined by phonons. Moreover, it is not 
clear what happens at the interface. Are the electronic and 
phononic systems in equilibrium, e.g. at the same tempera¬ 
ture, at the interface? If not are there any consequences for 
the interpretation of experimental results? 

To answer these questions, first we calculate the phonon 
thermal conductance k aitj and second we calculate the elec- 


TABLE I: Experimental thermal properties of bulk Fe and MgO. 


Thermal property at 300 K 

bcc-Fe fcc-MgO 

Phonon thermal conductivity k p (W/m/K) 

17.55 [20] 49.9 [21J 

Electron thermal conductivity k e (W/m/K) 

62.75 [20] 

Volumetric specific heat C (lO'T/irri/K) 

3.345 [20] 

Average phonon group velocity (m/s) 

7000 [22] 

Phonon Mean free path (nm) using Eq. 

6.39 
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FIG. 1: Sketch of a Fe/MgO(3 MLs)/Fe-MTJs used for phonon trans¬ 
mission calculation. The scattering region for the supercell calcula¬ 
tion is indicated. 
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tron thermal conductance k^ itj . Consequently, the overall 
thermal conductance k mtj is defined and the temperature 
profile across the Fe/MgO/Fe-MTJs is calculated by estimat¬ 
ing the electron-phonon coupling. 

The phonon transmission function and eventually the cor¬ 
responding contribution to the thermal conductance is calcu¬ 
lated by using the atomistic Green’s function (AGF) method, 
where details are described in Ref. [24] and references therein. 
In order to calculate the phonon transport the junction is di¬ 
vided into three parts shown in Fig. Q] The scattering re¬ 
gion consists of two Fe/MgO interfaces sandwiched by two 
semi-infinite Fe leads. After constructing the Green’s function 
G(e), the phonon transmission function t p {e) and the thermal 
conductance n p is calculated by [24] 


t p {e) = 


1 


(2 ?r ) 2 
Kp — — j £ t p (e) 


Tr[T L (e)G(e)T R (e)G\e)]d qil 

2DBZ /^\ 

df P {e,T )_ 


dT 


-de 


where T^e) and r K (s) are the broadening matrices of the 
bulk Fe leads, f p (e,T ) is the occupation function, and gy 
is the wave vector in the two dimensional Brillouin zone 
(2DBZ). The Green’s function is calculated from the inter¬ 
atomic force constants (IFCs). Thereby, we construct the 
Green’s function of the whole MTJ from the IFCs of bulk Fe 
and the IFCs of the scattering region. 

In order to obtain the IFCs we use density functional per¬ 
turbation theory (DFPT) li25ll implemented in the Abinit pack¬ 
age [26]. The self-consistent and perturbation calculations 
were performed by using Troullier-Martins norm-conserving 
pseudopotentials [27] and GTH LDA exchange-correlation 
functional [ 28[. First, to check the method we compute the 
phonon band structure for bulk Fe (a=2.867 A) and bulk MgO 
(a=4.239 A). For the self-consistent calculation of bcc Fe (fee 
MgO), we use an energy cutoff of 30 Ha (34 Ha) and a k- 
point mesh of 16 x 16 x 16 (8 x 8 x 8). For the dynamic 
matrix we use for both materials a 8 x 8 x 8 q mesh. The cal¬ 
culated phonon dispersions are shown in Fig. [2] In addition, 
experimental values are shown for comparison and we find an 
excellent agreement. 

To obtain the IFCs of the scattering region, a Fe/MgO/Fe 
supercell method is used. Thereby, the supercell consist 
of the MgO barrier with 4 monolayers of Fe at each side. 
The in-plane lattice constant is fixed to the value of bulk- 



FIG. 2: Phonon dispersions of (a) bcc-Fe and (b) fcc-MgO along high 
symmetry directions calculated by Abinit (color lines). The dotted 
data are experimental values taken from Refs. (2^1 and & Lines 
with different color correspond to different phonon branches. 


Fe (a=2.867 A). The cell volume and atomic coordinates are 
fully relaxed until the forces on all atoms are smaller than 
0.05 meV/A. After the structural relaxation, self-consistent 
and DFPT calculations are conducted according to different 
irreducible perturbations. For the self-consistent calculations 
we use a 10 x 10 x 2 k point mesh and an energy cutoff of 34 
Ha. For DFPT calculations, we use a 5 x 5 x 1 q mesh with 6 
inequivalent q points in the Brillouin zone. There are in total 
120, 150, 180 and 210 independent DFPT calculations for 3, 
5,7, and 9 monolayers of MgO, respectively. 

The phonon density of states (PDOS) projected on each 
atom in the scattering region is shown in Fig. da). The PDOS 
for the Fe layer near the MgO interface is significantly dif¬ 
ferent from the bulk PDOS of Fe due to the interface phonon 
bonding modes between Fe-MgO. This pronounced interface 
phonon peak is located at 8.5 meV. Inside the MgO barrier 
there are several phonon peaks. However, compared with the 
bulk phonon modes of MgO, the peak positions inside the 
MgO barrier shift to different energies. The phonon modes 
at the Fe-MgO interface and inside the MgO barrier may take 
effect through the electron-phonon interaction and they can be 
observed by inelastic electron tunneling spectra (IETS) [[31!]- 
For example, an IETS peak at 80 meV in d 2 I/dV 2 curve 
has been attributed to the Mg-O surface phonon modes in 
MgO|3.1] and this value is close to the largest peak position 
present in our projected PDOS. 

After obtaining the dynamic matrix, the IFCs in the junc¬ 
tion region are constructed. Away from the interface, the force 
constants are assumed to be that of bulk Fe. Our calculated 
phonon transmission function for different MgO thicknesses 
are shown in Fig.[3}b )- The cut-off energy of the transmission 
is at 38 meV and corresponds to the cut-off frequency of bulk- 
Fe. The general shape of the transmission function is almost 
independent of the MgO thickness. Only some shifts of the 
peak positions are visible, which may come from quantized 
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FIG. 3: (a) Projected phonon density of states (PDOS) for a supercell 
with 8 MLs Fe and 3 MLs MgO. The atom index is the same as 
in Fig. |T] The PDOS for bulk Fe and MgO is shown for reference, 
(b) The phonon transmission as a function of energy for Fe/MgO(3-9 
MLs)/Fe-MTJs.(c) The phonon and (d) electron thermal conductance 
of Fe/MgO(3-9 MLs)/Fe-MTJs as a function of temperature. 


phonon modes within the MgO barrier. Using the phonon 
transmission, the phonon thermal conductance of the MTJs 
Kp ITJ can be calculated by using eq. (0 and its temperature 
dependence is shown in Fig. 0c). There is a slight difference 
of Hp TJ between the different MgO thicknesses. At 300 K 
Kp ITJ is found to be at the order of 10 8 W/ m 2 / K. 

The screened-KKR Green’s function method 0 is 
adopted to calculate the electronic transmission t e (e) (details 
of the computation method are described in Refs. 113, 14]). 
k mtj j s 0 bt a i nec | f r0 m t e {e) by 132, 33 1 


-MTJ 1 rr (£l) 2 i 

-T [i2 ~“Lr ] 


(3) 


where L n = —— fj) n ^-f e (e)de, ji is the electro¬ 
chemical potential, and f e is the Fermi-Dirac distribution 
function. t e (e) depends on the relative magnetic orientation 
of the Fe leads to each other. 

The calculated k j ^ tj when the magnetization of the two Fe 
electrodes are in parallel or antiparallel alignment are shown 
in Fig. 0d). Since the electrons are tunneling through the 
MgO, k^ itj is quite small and decreases exponentially with 
increasing barrier thickness. Even at the smallest considered 
MgO thickness of 3 monolayers k:( v,7,/ at 300 K is 1 order of 
magnitude smaller than Hp ITJ . 

After obtaining Kp IT J and k ^ itj , the temperature profile 
in Fe/MgO/Fe-MTJs can be calculated by solving a one di¬ 
mensional thermal transport equation with appropriate bound¬ 
ary conditions. We consider the electron and phonon trans¬ 


TABLE II: Parameters used for calculating the electron-phonon in¬ 
teraction factor of Fe. 


A 

do (K) 

D{ef ) (States/Ha) 

G ep h (W/m 3 /K) 

0.243 (37] 

470 (38] 

23.02 (37] 

9.925 x 10 1 ' 


port, which carry the heat flux across the junction as well as 
the electron-phonon interaction in the Fe leads. It is reason¬ 
able to neglect the magnon contribution to the thermal trans¬ 
port since the magnon transport across the junction is assisted 
by the electron tunneling through electron-magnon interac¬ 
tion, which is at least one order of magnitude smaller than the 
electron and phonon transport process. The two fluid model 
which take into account the electron-phonon non-equilibrium 
in the Fe electrodes will result in different electron and phonon 
temperature T e (z) and T p (z) at the Fe-MgO interfaces. 

The energy balance equations between electrons and 
phonons in the Fe leads are 0 

k FefW_ _ Geph (T e -Tp) = 0 

( 4 ) 

k p e ^£^- + G eph(Te -T p ) = 0 

where /cf e and k p e are the electron and phonon thermal con¬ 
ductivities in Fe, respectively. For these quantities we take 
the bulk values, which are listed in Table U G ep h is the 
electron-phonon interaction factor, which is determined by 
G ep h = TrkBM u}2 )D( £ F) [35]. In this expression, kp is the 
Boltzmann constant, A is the electron-phonon mass enhance¬ 
ment parameter, D{ep) is the electron density of states at the 
Fermi energy, and ( oj 2 ) is the second moment of the phonon 
spectrum defined by McMillan 0. (co 2 ) can be approxi¬ 
mated by (w 2 ) ps 9 2 d /2 where dp is the Debye temperature of 
Fe. By using the related parameters listed in Table. ITTlGAoj, of 
Fe is found to be 9.925 x 10 17 W/m 3 / K. This value for Fe 
is comparable but slightly smaller than that of Nickel (10.5 x 
10 17 W/m 3 /K) and Platinum (10.9 x l0 17 W/m 3 /K) 0. 

The solutions of Eq. 0 for the left (L) and right (R) Fe 
layers are 



where the coefficients of the electrons and phonons hold 


B L(R) = B L(R) . C L(R) 


h.Fe 

C L(R) . n L(R) _ _'^P_ D L(R) 
* e k Fe P 


Thus the last term in Eq. 0 accounts for the electron-phonon 
imbalance at the Fe/MgO interface. Consequently, 


iFe = 


1 uFeh,Fe 

1 


G eph fcf 6 




( 6 ) 


is the characteristic length of this imbalance. We define 


sFe 


G, 


eph 


l/lle = Vll + Vll\ 


(7) 
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In total we have to calculate 6 coefficients B r J' R> , Cl JR \ and 
L( R) 

D e y by using the following 6 boundary conditions 


dT e,ph( Z ) , 

dz U= - d/2 

MTJ ( t L (Z^L\^T R <-V\ 
^e^ph 1 1 e,ph\ ^ e,ph\ J 

T e L (z = -L) = T l 


dT e R ph (z) 

dz 


z—dj"! 


Fe 


dT e R ph (z) 


n ‘e,ph \z=d/2 

T r (z = L) = T r 


where the last two conditions are some given temperatures at 
left and right of the MTJ, which has a total thickness of 2 L. 

To illustrate our result we show in Fig. Q] the electron 
and phonon temperature profiles for a Fe(10 nm)/MgO(9 
MLs)/Fe(10 nm)-MTJ with a temperature difference of 1 K. 
The consequence of different electron and phonon tempera¬ 
tures near the MgO/Fe interface is significant especially for 
the definition and evaluation of thermoelectric physical quan¬ 
tities. For example, the Seebeck coefficient S is defined as 
S = where A T e is the electron temperature drop since 
the Seebeck effect is a consequence of electron transport. 

Although we have a non-equilibrium situation it might 
be useful to define an effective temperature drop AT e //, 
which can be used to define a k mtj = -rS —, where q = 

—k Re dT j i ^ ~ k p e dlp ,] { - ] is the total thermal current through 
the junction. Note that q is conserved whereas the electron and 
phonon thermal current alone are not due to the imbalance. 

Following Ref. [34] we define the effective temperature 
T e ff by linear extrapolation of the electron or phonon tem¬ 
perature from the equilibrium towards the MgO barrier (see 
Fig. [4]). For the case shown in Fig. [4] we calculate k aitj = 
1.681 x 10 8 lF/m 2 /iv. This value can be used in a simple 
network model to estimate the electron temperature drop by 
calculating A T e ff. To get an idea about the error in this esti¬ 
mation we calculate the following ratio 

Te T eff 

Teff - Tp k?‘ 

This means if k Re k Re , which is the case for most metals, 
T e ff is closer to the electron temperature than to the phonon 
temperature. If this is not the case or a more precise knowl¬ 
edge of A T e is needed one has to solve appropriate transport 
equations, i.e. Eq. ©, for the whole junction, but taking 
the first principle values for k^ [tj and ti p TJ . Moreover, the 
effect of the imbalance is larger if the phonon interface con¬ 
ductance is larger. This effect is shown in the inset of Fig. [4] 
where we assume a ten times larger n p ITJ than the first princi¬ 
ple value, which can be achieved by using different materials. 
Note that even the phonon temperature drop across the barrier 
decreases the drop in the electron temperature remains. 

The thermal conductance n MTJ of a MTJ is of particular 
importance for the interpretation of experimental results. In 
experiments a Seebeck voltage is measured and a temperature 
drop is estimated using diffusion models. As stated earlier 
often a value of bulk or thin film value of MgO is used for 


-L -d/2 d/2 L 



FIG. 4: The electron temperature T e , the phonon temperature T p and 
the linear extrapolation of electron temperature Te// profiles across a 
Fe/MgO(9 MLs)/Fe-MTJ with the temperature at the left side (L=-10 
nm) Tl = 30Ilf and right side (L=10 nm) Tr = 30010 The inset 
shows the temperature profiles assuming a ten times larger n p TJ 
than the first principle value. 



d (nm) 

FIG. 5: The junction thermal conductance as a function of MgO 
barrier thickness in this work (black) compared to bulk MgO (red 
dashed) and thin film MgO (blue dashed). The value used by Walter 
et.al. 01 is shown as a green star. 


k mtj . In Fig. [5]we plot our calculated n MTJ in compari¬ 
son to bulk and thin film values of MgO for different MgO 
thicknesses. Our results are almost independent of the MgO 
thickness. Further, as an example the value used in Ref. [10] 
is indicated by the green star. This implies that the estimated 
Seebeck coefficients in Ref. El may be too high. 

In conclusion, we show that there is an imbalance of the 
electron and phonon temperature at nano-magnetic interfaces. 
This leads to different temperature drops for electrons and 
phonons and one has to be cautious interpreting results. Nev¬ 
ertheless, in the considered Fe/MgO/Fe MTJ we observed a 
large interface resistance for the phonons. In consequence, 
a large drop also for the electron temperature exists, which 
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is responsible for high Seebeck voltages observed in experi¬ 
ments. But even for MTJs with low phonon interface resis¬ 
tances the drop of the electron temperature remains large due 
to the strong imbalance of electron and phonon temperature. 
For a MTJ with metallic leads we define an effective thermal 
conductance by defining an effective temperature. The corre¬ 
sponding values given in Fig. [3 can be used in simple mod¬ 
els as the conductance of the MTJ (barrier plus interfaces). 
However, if more reliable values of the temperature drops are 
needed a calculation of the coupled electron and phonon sys¬ 
tem is necessary. In any case for the qualitative description 
first principle calculations are a need due to the coherent trans¬ 
port across the nano-magnetic interface. 
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